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Non dispersive electronic Rydberg wave packets may be 
created in atoms illuminated by a microwave field of circu- 
lar polarization. We discuss the spontaneous emission from 
such states and show that the elastic incoherent component 
(occuring at the frequency of the driving field) dominates the 
spectrum in the semiclassical limit, contrary to earlier predic- 
tions. We calculate the frequencies of single photon emissions 
and the associated rates in the "harmonic approximation" , i.e. 
when the wave packet has approximately a Gaussian shape. 
The results agree well with exact quantum mechanical calcu- 
lations, which validates the analytical approach. 

PACS:32.80.-t, 42.50-p, 32.50.+d, 03.65.Sq 



I. INTRODUCTION 

Recent theoretical studies predicted the existence of 
non-dispersive wave packets in atoms driven by a mi- 
crowave radiation of either circular ]IJ or linear polariza- 
tion ||. Contrary to standard electronic wave packets 
which disperse in time when the binding force is nonlin- 
ear, these novel objects do not spread due to a nonlinear 
coupling between the atom and the driving field which 
locks the electronic motion onto the driving frequency. 
The phenomenon responsible for this locking is a 1:1 res- 
onance between the internal (Kepler) and external (mi- 
crowave) frequencies, creating a resonance island at the 
center of which lies a stable periodic orbit. In the dressed 
atom picture || , a single eigenstate of the full "dressed" 
hamiltonian represents the wave packet [^]j4j, explaining 
immediately the absence of spreading. 

The very same microwave field is also responsible for 
the decay of the wave packet via ionization. As shown 
elsewhere ||, the typical ionization time of wave pack- 
ets built from circular Rydberg states around principal 
quantum number ng — 60 exceeds 10 5 Kepler periods be- 
ing often longer than 10 6 periods. The ionization occurs 
via a "chaos assisted tunnelling" mechanism || (the tun- 
neling occurring from the stable island into the chaotic 
sea surrounding it). In effect the lifetime against ioniza- 
tion is very sensitive to small changes of the parameters 
(microwave frequency, f2 or its amplitude F) exhibiting 
large scale fluctuations ||. 

Another possible decay mechanism of the wave packet 
is of course the spontaneous emission of photons. For 
circularly polarized microwave, the situation considered 



in this paper, it has been discussed in Here we read- 
dress this issue since, in our opinion, the results reported 
in JjJ are incomplete. To this end, we use a different 
approach, borrowing heavily from the well known treat- 
ment of resonance fluorescence of two-level atoms driven 
by a laser field jj| using the dressed atom approach. The 
main advantage of the present approach is its generality: 
it can be used for any polarization of the microwave field, 
for any superposition of driving fields with various fre- 
quencies and for any structure of electromagnetic modes 
surrounding the atom (as shown e.g. in the presence 
of a cavity can strongly affect the spontaneous emission 
of an atom) . Also, it uses the common knowledge widely 
developed in the quantum optics community and is thus 
accessible to a large audience. Let us mention also that 
spontaneous emission of non-dispersive wave packets in 
linearly polarized microwave is discussed in Q . 

The paper is organized as follows. In section II, we 
discuss the basic properties of the non-dispersive wave 
packets (not considering spontaneous emission) from dif- 
ferent points of view: laboratory frame or rotating frame, 
time-dependent or Floquet pictures, dressed atom pic- 
ture. The discussion is necessary for a proper under- 
standing of spontaneous emission, which is done in sec- 
tion III. Our results are compared to the earlier published 
results in section IV. Section V is devoted to the valid- 
ity of the harmonic approximation, by comparing our 
analytical results to "exact" numerical calculations. We 
conclude in section VI. 



II. THE NON-DISPERSIVE WAVE PACKETS 
FROM DIFFERENT POINTS OF VIEW 

We consider a hydrogen atom in the presence of a 
monochromatic driving microwave field of circular po- 
larization, that is an electromagnetic field which can be 
described by a single mode. In this section, we will not 
consider any coupling to other (initially empty) modes of 
the electromagnetic field. This will be discussed in sec- 
tion III. We first treat the microwave field classically. A 
full quantum tr eatm ent of the microwave field is possi- 
ble (see section 



IIP ), but not really needed as, for usual 



microwave fields, the number of photons in the mode is 
enormous and a classical treatment of the microwave field 
is possible. 
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A. Hamiltonian in the laboratory frame 

With the previous assumption, the quantum Hamil- 
tonian of the system is (in atomic units, neglecting rel- 
ativistic and QED effects, assuming an infinite mass of 
the nucleus and in the dipole approximation): 

p 2 1 

H c i = h Flxcosttt + ysiaClt). (1) 

2 r 

where is the frequency of the microwave field circularly 
polarized in the x — y plane and F its amplitude. The 
x-axis is chosen to coincide with the microwave electric 
field at t = 0. 



B. Floquet approach 

The Hamiltonian ([!]) is periodic in time with period 
27r/0. Hence, a Floquet approach |l(J may be used: any 
solution of the time-dependent Schrodinger equation can 
be expanded as a linear combination of time-periodic 
states - called Floquet states of the system - with co- 
efficients oscillating in time as exp(— iEtf) where Ei is 
called the quasi-energy of the Floquet state (there is no 
explicit h in the equations as we are using atomic units 
where it is one). Obviously, the quasi-energy is defined 
only modulo Q making the Floquet spectrum periodic in 
energy with period precisely f2. The Floquet spectrum as 
well as the Floquet states can be obtained from the di- 
agonalization of a suitable operator, namely the Floquet 
Hamiltonian: 

H = H d - (2) 

with time-periodic boundary conditions. A suitable basis 
for the Floquet Hamiltonian is the product of any basis in 
usual space - for example, the usual hydrogenic eigenbasis 
with bound states \nlm) - by the basis of time-periodic 
oscillating exponentials exjp(iKflt) with K an integer. 
We will denote \nlm; K) such a basis. The next step is 
to expand the time-periodic Hamiltonian in a Fourier se- 
ries and calculate the non-zero matrix elements of the 
Floquet Hamiltonian. Because no higher harmonics of 
the microwave frequency exist in the Hamiltonian, only 
matrix elements changing K by zero or one unit do not 
vanish. Moreover, in the specific case of circular polariza- 
tion, the part in the Hamiltonian which increases (resp. 
decreases) K by one is proportional to x—iy (resp. x+iy) 
hence decreases (resp. increases) the angular momentum 
of the electron along the z axis by one unit. In other 
words, the Floquet Hamiltonian obeys the selection rule 
A(K + m) — 0. This signifies merely that absorption of a 
o~ + polarized microwave photon must increase the atomic 
angular momentum L z — xp y — yp x by one unit. The Flo- 
quet Hamiltonian thus has a block diagonal structure and 
the Hilbert space can be split into manifolds character- 
ized by constant value of K. = m + K. In each block, the 



Floquet Hamiltonian is no longer periodic in energy. On 
the other hand, the Floquet Hamiltonian in blocks JC and 
KJ only differ by a constant energy term (JC' — K)Q which 
restores the total periodicity of the Floquet spectrum. 
However, diagonalization in a given block is sufficient to 
obtain all the information on the system. 

C. Hamiltonian in the rotating frame 

For a circularly polarized radiation, another possibil- 
ity exists which has been frequently used in the past 
By passing to the non inertial frame rotat- 
ing with frequency Q using the unitary transformation 
U = cxp(iQL z t), one may remove the oscillatory time 
dependence from Eq. (|l]) obtaining in the rotating frame 
the Hamiltonian 

dW p 2 1 
H Iot = U H cl f/t + iU— = ^---+Fx- QL Z . (3) 
ot 2 r 

This time-independent Hamiltonian has some energy 
levels and corresponding eigenstates. Obviously, its spec- 
trum is not f2-periodic. However, the link with the Flo- 
quet approach in the laboratory frame is easily estab- 
lished: there is a one-to-one correspondence between the 
eigenstates in the rotating frame and the Floquet states 
in the laboratory frame from the block fC — 0. Indeed, 
consider an eigenstate (i.e. time-independent state) \4>i) 
of H Iot with energy Ei. From Eq. we deduce straight- 
forwardly: 

\H cl -i^\u^\^) = E i U^^) (4) 

which implies (see Eq. (|J)) that U^\<pi) is a Floquet eigen- 
state with quasi-energy Ei. From the definition U' = 
exp(— iflL z t), it is clear that the K-th Fourier compo- 
nent of U' \(f>i) is nothing but the component of \<j>i) on 
the subspace with L z = m = —K. In other words, £^|</>;) 
is a Floquet eigenstate in the laboratory frame lying in 
the block K. = m + K = 0. 

Of course, Floquet states belonging to the other blocks 
K / can also be found from the time-independent 
Hamiltonian in rotating frame. A first method is to use, 
instead of the operator U, a slightly different unitary 
equivalent operator exp(— ikflt)U (with k an integer, it 
has the same time-periodicity as U itself) which leads to 
the following Hamiltonian in the rotating frame: 

H Iot = 2- - - + Fx - Q(L Z - k). (5) 
2 r 

Taking all possible positive and negative integer k values 
one recovers the full Floquet spectrum. 

A second equivalent method is to keep the standard 
Hamiltonian in rotating frame, Eq. (^|), and to consider 
states of the form exp(ikttt)W \<fii) . These are time- 
periodic Floquet (if k is an integer), eigenstates of the 
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Floquet Hamiltonian H., Eq. (g) with energy Ei + kVL, 
belonging to the block K = m + K = k. 

At a deeper level, this can be seen as a consequence 
of the unitary equivalence (which implies the same spec- 
trum) for the Floquet Hamiltonians in the laboratory 
frame and in the rotating frame. The spectrum of the lat- 
ter is here very simple: because H ro t is time- independent, 
it is simply the energy spectrum of H mt shifted by an ar- 
bitrary integer multiple of ft. For a more complicated 
case like - for example - a microwave with polarization 
close to circular, the Hamiltonian in the rotating frame 
is weakly time-dependent and its corresponding Floquet 
Hamiltonian (with the time dependent term treated per- 
turbatively) may be the simplest approach. 



D. Dressed atom approach 

An alternative view of this problem is to treat the mi- 
crowave field quantum-mechanically as being in either a 
coherent state or a pure Fock state with a large num- 
ber of photons in the circularly polarized mode p2[ . 
This approach is quite common in quantum optics and 
is known as the dressed atom picture. Then the atom is 
not considered as an isolated quantum system, but it is 
rather the strongly coupled atom+field system which is 
treated using quantum mechanics. The Hamiltonian of 
the atom+field system is: 



+ ttah 



(6) 



where a^/a are the creation/annihilation operators in the 
occupied, a + circularly polarized mode, d± = (xiiy) / y/2 
the corresponding components of the dipole moment op- 
erator and V is the coupl ing st rength, whose exact ex- 
pression is given in section III A but not important at the 
present level. 

The Hamiltonian H q is time independent and may be 
expanded in the basis: {\nlm) \N) — \nlm;N}} where 
\nlm) are hydrogenic states and N denotes the number of 
microwave photons in the occupied mode. The diagonal- 
ization of the Hamiltonian yields then the dressed states; 
as its stands, the energy spectrum is not fi-periodic. 

The absorption of a circularly polarized microwave 
photon increases by one the projection L z of the elec- 
tron angular momentum on the z axis. In the quantum 
field basis this manifests itself in the fact that the Hamil- 
tonian matrix is block diagonal. The Hilbert space splits 
into manifolds characterized by the total atom+field pro- 
jection of the angular momentum A/", i.e., basis states 
coupled within a given block have fixed Af = m + N. Af 
is simply the total angular momentum projection on the 
z axis carried by both the electron and the photons. 

When the number of photons is very large, N 3> 1, we 
can neglect the variations in the non-diagonal matrix el- 
ements - proportional to either y/N or y/N + 1 - across 
neighboring N blocks. Then, >/2Vy/N can be identified 



as the amplitude of the microwave field F. With this ap- 
proximation, the energy spectrum is (locally) Sl-periodic. 
As a matter of fact, the matrix elements of H q are ob- 
viously equal to the ones of the classical Floqu et H amil- 
tonian, Eq. (||), already discussed in section [IB, pro- 
vided that the quantum number K labelling the Fourier 
components in the Floquet picture is identified with the 
number of photons N. This establishes - in the limit of a 
large number of photons - the one-to-one correspondence 
between the quasi-energy spectrum of the Floquet Hamil- 
tonian and the energy spectrum of the dressed atom, and 
also the correspondence between the Floquet eigenstates 
and the dressed atom eigenstates. The two points of view 
are completely equivalent: in the following, we will use 
the dressed atom picture as it seems more convenient for 
incorporating the spontaneous emission of photons. 

As a final remark, let us note that the previous dis- 
cussion is already well known in a different context in 
quantum optics: for the simple case of a two level atom 
coupled resonantly by an oscillatory field, in the rotating 
wave approximation. There, a field-quantized or Floquet 
approach yields an energy spectrum periodic in f2 - a se- 
quence of doublets separated by the Rabi frequency || 
- while semiclassically in the "rotating frame" (2^] one 
gets a single doublet of states. 



E. Numerical calculation of the energy spectrum 

As discussed in the three preceding sections, there ex- 
ists a one-to-one correspondence between the eigenstates 
in the rotating frame, the Floquet eigenstates and the 
dressed states in the quantum language. Therefore, we 
shall call in the following the cigcncncrgies of either of the 
Hamiltonians (||),(|3|) or (^) as quasienergies while refer- 
ring to the eigenstates as Floquet eigenstates or dressed 
states. 

As we are interested in the regime of highly excited 
Rydberg states (with typical principal quantum number 
60) driven by a resonant microwave field, the number of 
atomic states significantly coupled by the microwave field 
is enormous. For computational efficiency, it is better to 
use a Sturmian basis |J-|6| rather than the usual atomic 
basis. The advantage is that it provides a discretization 
of the atomic continuum. Combined with the complex 
rotation technique, it allows us to find not only the ener- 
gies of the dressed states but also their lifetimes against 
ionization by the microwave field ]5|,[|. To numerically 
find several converged dressed states, diagonalizations of 
very large matrices (of the order of 10 5 ) are necessary. As 
discussed above, there exists a full equivalence between 
diagonalization of the Floquet Hamiltonian, i.e. clever 
use of the time-periodic structure of the problem and di- 
agonalization of the Hamiltonian in the rotating frame. 
The computational efficiencies of the two points of view 
are strictly equal. Our computer code ||-p|_diagonalizes 
in fact H TOt , contrary to the statement in Eg]. 
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F. Non-dispersive wave packets 

Some time ago, Klar [^4| noticed that the Hamiltonian 
(|J) in the rotating frame allows for the existence of a 
stable fixed (equilibrium) point in a certain range of mi- 
crowave amplitude F . The beautiful contribution of |l|] 
was to notice that wave packets initially localized in the 
vicinity of this fixed point will not disperse (being bound 
by the fact that the fixed point is stable) for very long 
time. In the laboratory frame, these wave packets also 
called "Trojan states" in JlJ appear as wave packets mov- 
ing around the nucleus along the circular trajectory [the 
periodic orbit of (|l|)]. The existence of "stationary" wave 
packets - keeping exactly their shape during propagation 
along the circular orbit and hence not dispersing at arbi- 
trarily long times (neglecting the very slow ionization) - 
was later found in [|J . These are simply some among the 
eigenstates of the dressed atom: by construction, they 
are strictly periodic functions of time. 

Semiclassically, the wave packets are manifestations of 
the motion inside the stable primary 1:1 resonance zone 
when the microwave frequency matches the electronic 
frequency. Similar states have been studied earlier for 
model one-dimensional driven systems [ p5| . While for 
other models, the wave packet resembles itself every pe- 
riod of the driving perturbation (while changing shape 
within one period ^,^5|), the wave packets in circular 
polarization keep the same shape at all times. This is 
due to the conservation of Af. 

Quantum mechanically, any dressed eigenstate may be 
viewed as a linear combination of atomic states, com- 
bined with appropriate photonic states (in such a way 
that in the combination Af is conserved). The wave 
packet states are linear combinations of atomic Ryd- 
berg states coupled by resonant microwaves. Thus, the 
biggest contribution to this combination comes from cir- 
cular states with principal quantum number n close to: 

n = rr 1 / 3 . (7) 

Note that no is not necessarily an integer. For the partic- 
ular case of the "ground state wave packet" (with maxi- 
mum localization near the fixed point), it is a superposi- 
tion of mainly circular atomic states \n, n — 1, n — 1) with 
a small contamination of other quasi-two-dimensional 
states of the form \n,l,l) and even smaller contributions 
from states I ^ m extending significantly out of the xy 
plane. 

The fixed point in the rotating frame is located at y e = 
z e = and x e satisfying 

\/x\ + F - x e n 2 = 0. (8) 

It is convenient to introduce [Q a dimensionless parame- 
ter: 

q = fr 2 x~ 3 . (9) 



Then q = 1 for F — and the fixed point is stable 
for q e [8/9,1]. It is also useful to define scaled vari- 
ables which make the discussion of the results less depen- 
dent on the particular microwave frequency value [B6| . 
We have already defined the effective principal quantum 
number no, Eq. (0). The scaled microwave amplitude is 
defined as 

F = Fnt = MT 4/3 . (10) 

Fq is a dimensionless parameter simply related to q, see 
Eq. (|), via 

F = (1 - q)q- 1 '\ (11) 

A possible way for creation of a non-dispersive wave 
packet consists of the preparation of a circular state 
with a given principal quantum number no followed 
by a smooth turn on of the microwave with frequency 
n = 1 /rig • Such an approach has been tested numer- 
ically by integration of the time-dependent Schrodinger 
equation ]27|] . It has been shown that the wave packet 
may be obtained with about 90% efficiency. 

During the switch-on of the microwave, two regimes 
may be identified. For q only slightly less than unity, 
(small microwave amplitude), the resonant coupling be- 
tween circular states is important only. The effective ap- 
proximate Hamiltonian in the vicinity of the fixed point 
corresponds then to the Hamiltonian of a pendulum [fl2[ . 
This is the standard description of the resonance at first 
order (28). For low Fo, the resonance island is small in 
action variable and, in effect, the pendulum states must 
be quite extended in the conjugate angle. In the rotat- 
ing frame, this manifests itself as the fixed point being 
marginally stable along the y direction at vanishing Fq. 
For too small Fo, the binding force is too small and the 
wave packet eigenstate can be hardly considered to be 
localized in the very vicinity of the fixed point and, in 
turn, cannot be considered as a localized wave packet in 
the laboratory frame. 

With increasing microwave amplitude, the resonance 
island grows and may support a few (depending on the 
size of the effective fr = 1/no) states. Then, for a suffi- 
ciently large Fo (but not large enough to destabilize the 
resonance), the low excitations in the resonance island 
may be well described in the harmonic approximation. 
It is this regime of Fo values which is of interest. 

The harmonic expansion of the Hamiltonian iJ ro t in 
rotating frame around the equilibrium point allows one 
to find normal modes of the problem and obtain ana- 
lytic predictions concerning the energies of wave packet 
eigenstates |l],||. 

G. Harmonic approximation 

The Hamiltonian H Iot can be expanded in powers of 
position and momentum operators in the vicinity of the 
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fixed point. The first order terms cancel because this 
is a fixed point. In the harmonic approximation, only 
second order terms are kept (including of course crossed 
position-momentum terms) and we obtain: 



'H h =E e + 

r>2 ~2 
— ii qx 



P 2 r,/~ - \ Q 2 qy 2 
— + \l[xp y - yp x ) + — — 



Q 2 qP 



(12) 



where 



E - i^si 2 / 3 

e - 2g 2 /3 



(13) 

is the energy of the equilibrium point and [x, y, z) = (x — 
x e , y — y e , z — z e ) denotes the displacement with respect 
to the fixed point. 

The "diagonalization" of this Hamiltonian leads to the 
definition of the normal modes for the electronic motion. 
The calculation is tedious but straightforward: the fre- 
quencies of the 3 normal modes labelled by the (+, — , z) 
indices are: 



2-q±< 



with a shorthand notation 



Q = VV - 



(14) 
(15) 

(16) 



The corresponding creation and annihilation operators 
in the 3 normal modes are linear combinations of po- 
sition and momentum operators. In the harmonic ap- 
proximation, the motion perpendicular to the plane of 
microwave polarization decouples from the motion in the 
plane; hence, the annihilation operator in the z-mode has 
the following simple expression: 



q 1 / 4 VttS + iq 1 / 4 p z /Vn 



(17) 



The creation operator b\, being the adjoint of b z , has a 
similar expression with imaginary part of opposite sign. 
The ± modes completely entangle the various coordi- 
nates and momenta (x,y,p x ,p y ). They are neither the 
usual "cartesian" normal modes along x and y, nor "cir- 
cular" normal modes along x ± iy, but intermediate com- 
binations. The annihilation operators are: 



b ± =N ± (y/Ux-i {Q± Y m 



3/2 



AQ±q)VQ - , 3q±Q - 
2quj ± fx I" 2qy/U Py 



with the normalization constants given by: 
V(2Q 2 -g + 2±3(l-2<?)Q) 



N± 



2 (l-9) 



1/4 



(18) 



(19) 



The Hamiltonian in the harmonic approximation can be 
expressed in terms of creation and annihilation operators 
of the electronic normal modes [flljilJ 



(20) 



H h = E e +cj+(b\b + + l/2) 

+ 1/2) + uj z (blb z + 1/2). 



Note the unusual minus sign in front of the u>_ term. 
This comes from the fact that the stable fixed point is not 
a minimum of the Hamiltonian. However, the presence 
of crossed momentum-position terms assures the com- 
plete stability of the equilibrium point versus any small 
perturbation, in complete analogy with a usual potential 
minimum. 

Let us denote the atomic eigenstates in the harmonic 
approximation as \n +1 n-,n z ) where n, denotes the ex- 
citation in a given ±, z mode We will also - when 
needed - denote the Floquet or dressed atom eigenstates 
(thus including the microwave field degree of freedom) as 
n_, n z \ AT). The corresponding energies are: 

E(n + ,n-,n z ;N) = E e +N£l + (n+ + - J uj+ 



L 












- ( n z - 













It has been verified that exact Floquet quasienergies, 
corresponding to dressed states localized in the vicin- 
ity of the equilibrium point in the rotating frame agree 
with the harmonic prediction, Eq. (21), to within 10% 
of the mean level spacing |^,|| for wave packets corre- 
sponding to no w 60 or higher. This indicates that the 
harmonic approximation is quite satisfactory for suffi- 
ciently large uq (sufficiently small si). In particular, the 
"ground state", |0, 0,0) corresponds to the most tightly 
bound state, centered at the fixed point in the rotating 
frame. Its wavefunction has a Gaussian shape in the har- 
monic approximation. The corresponding exact Floquet 
(dressed) state || is a non-dispersive wave packet since in 
the laboratory frame it moves in time around the nucleus 
with frequency Q following the stable circular periodic 
orbit. It is precisely the spontaneous emission from this 
state which will be discussed in the following sections. 

The partial scheme of levels resulting from Hh is shown 
in Fig. |l| for a few values of the "total" angular momen- 
tum Af. The neighboring ladders of states are shifted 
with respect to each other by the microwave frequency 
il corresponding to different decoupled blocks of the Flo- 
quet Hamiltonian. Transitions between quasienergy lev- 
els (Floquet or dressed states) may be introduced by 
some additional perturbation. In particular one may con- 
sider the so-called Floquet spectroscopy Q , absorption 
or stimulated emission between dressed states introduced 
by an additional laser or microwave field. 

It seems quite natural to study the possible effects of 
spontaneous emission of these states within the very same 
harmonic approximation. This allows to obtain close an- 
alytic expressions for the spontaneous decay rate. The 
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validity of the harmonic approximation has been chal- 
lenged strongly by Farrelly and coworkers [^6|-^o|. By 
comparison of spontaneous spectra obtained in the har- 
monic approximation with full quantum mechanical re- 
sults, we show later that their claims are not justified. 

III. SPONTANEOUS EMISSION OF THE WAVE 
PACKET IN THE DRESSED ATOM PICTURE 

A. General properties 

To consider the spontaneous emission one has to con- 
sider the atom interacting not only with a single mode 
[treated either classically (|l|) or quantum mechanically 
@] but a full electromagnetic field. The corresponding 
Hamiltonian reads, in atomic units |30|j 

H = Hi + r • E(0) + H[ (22) 

where subscript i stands for cl or q corresponding to the 
classical or quantum version of the description of the 
occupied microwave mode, respectively, r is the usual 
dipole operator, E(r) denotes the electric field opera- 
tor at position r and Hf is the free field Hamiltonian. 
Both Hf and E can be expanded over modes of the 
field. While in §, the expansion in sperical waves is 
used (which facilitates explicit transformation to the ro- 
tating frame), we expand the field into standard plane 
waves H]. We normalize the modes per unit in k space 
(related to the frequency in atomic units by lu = ak 
where a — e 2 /Tic ~ 1/137.036 is the fine structure con- 
stant), then 

E(r) = ^ f ^-(a j (k)e j exp{ik-r)+h.c.) d 3 k, (23) 

J 67T 

3 

where h.c. stands for hermitian conjugation, Oj(k) is 
the annihilation operator in the mode characterized by 
wavevector k and polarization &j (the sum over j is over 
the two possible orthogonal polarizations). Note that 
in ( p2| ) the electric field operator is evaluated at origin, 
consistently with the dipole approximation assumed in 
this paper. Similarly the free field hamiltonian is 

H f = ^ ( w «j( k ) ta j( k ) d3k - ( 24 ) 
j J 

Now we can use the Fermi golden rule to calculate the 
decay rate due to spontaneous emission from an arbitrary 
initial state \i) of the atom dressed by the microwave field 
to any final state |/), taking into account the density of 
modes for the photons fl3l| . The result is: 

r=^|(/|r|z)| 2 (25) 

where lo is the energy difference between the initial and 
final dressed states. 



Spontaneous emission will take place between the var- 
ious states belonging to ladders depicted in Fig. [l| In 
the process, the projection of the full angular momen- 
tum on the Oz axis is conserved. Here, this total angular 
momentum has three contributions: the atomic angular 
momentum m, the angular momentum N of the driving 
microwave photons (all in the same mode) and the angu- 
lar momentum M in the other modes of the electromag- 
netic field. Hence, we have A4 — Af+M = N + m + M = 
constant. As spontaneous emission is here considered as a 
weak perturbation, it will simply induce transitions be- 
tween eigenstates of the dressed atom, see Eq. (^|), for 
which TV is a good quantum number (see preceding sec- 
tion). In the dipole approximation, for a one-photon pro- 
cess, the atomic angular momentum m can change by at 
most one. When the spontaneous photon is emitted in an 
initially empty mode (which is the process we are inter- 
ested in) , the number of photons in the driving microwave 
mode N does not change. Hence, TV changes by at most 
one, and the same is true (with opposite sign) for M. 

Thus, in the dipole approximation, there are only three 
possible processes. If the spontaneous photon is emitted 
with <j + polarization, i.e. angular momentum M = 1, it 
will induce a transition from a state in a given TV ladder 
to a state in the Af' — Af — 1 ladder. Similarly, for a~ 
polarization, i.e. angular momentum M = —1, it will 
induce a transition from Af to TV' = TV + 1 and for it 
polarization, M — and TV' = TV. As all but one photon 
modes are initially empty, by spontaneous emission, the 
energy Two is gained in one mode and thus has to be lost 
by the dressed atom. Note that here and in the rest 
of this paper, we discuss the polarization of the various 
transitions using a quantization axis along the z axis, 
the direction of the propagation of the driving field. As 
explained below the spontaneous photons are mainly a + 
polarized with respect to the z axis. Hence, the helicity of 
the photons (with respect to their direction of emission, 
i.e. the quantity measured in an experiment) is in general 
elliptical, depending on the direction of emission. 

In summary, in the dressed atom picture as shown in 
Fig. [j], the only possibilities for emitting a spontaneous 
photon is to decay to a state with lower energy in the 
same ladder or in one of the two neighboring ladders. All 
other processes are forbidden. Note however that - in the 
absence of a further approximation - plenty of processes 
with all 3 possible polarizations should exist. However, 
as can be seen from Fig. [I], the cr + emission is clearly 
favored, since the TV— 1 ladder is shifted downwards with 
respect to the initial state by Q. n is less favored and a~ 
even less because of the upwards energy shift. 
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FIG. 1. Schematic energy levels of an hydrogen atom 
dressed by a microwave field of circular polarization, with 
Hamiltonian given by Eq. (jf]). The total angular momentum 
Af (atom+field) is a constant of motion, so that the spectrum 
can be split in separate ladders. The consecutive ladders have 
the same structure, simply separated (in energy) by the mi- 
crowave frequency Q. We have here represented three consec- 
utive ladders. If all ladders are taken into account, the spec- 
trum is f2-periodic and coincides with the Floquet spectrum 
of the classical Hamiltonian, Eq. (^). Alternatively, the spec- 
trum in each ladder can be thought as the energy spectrum of 
the Hamiltonian in rotating frame, Eq.Q. Spontaneous emis- 
sion of photons (represented by arrows in the figure) from a 
state in ladder Af is only possible to another state lying at a 
lower energy in either ladder Af — 1 (with a + polarization) , 
ladder Af (with n polarization) or ladder Af + 1 (with cr_ po- 
larization). Hence, a + transitions are favored. In the specific 
case of the ground non-dispersive wave packet |0, 0, 0) in the 
harmonic approximation, only three transitions, represented 
by the solid arrows, are allowed at frequencies fi — u+ , il and 
Q + lu- , all of them being a + polarized. 



B. Spontaneous emission from a non-dispersive wave 
packet in the harmonic approximation 

The general properties discussed in the preceding sec- 
tion are valid whatever initial state is chosen. When the 
initial state is a non-dispersive wave packet and when the 
scaled microwave field F is sufficiently large, we can use 
the harmonic approximation described above. We will 
make this approximation throughout the rest of this sec- 
tion. The atomic dipole components x^iy,z correspond- 
ing respectively to the photon polarizations a± , ir can be 



expressed as functions of creation/ annihilation operators 
in the normal modes, Eqs. (^7|)-(^8|): this is straightfor- 
ward as all these operators are linear combinations of 
positions and momenta. The result is: 



x ±iy 




-3q + Q±(q-Q) 
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3q + Q T {q + 
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+ 

<T 1/4 + 





(26) 



(27) 



Hence, in the harmonic approximation, there are only 
very few transitions allowed. The selection rules are: 



An+ = 
An + = 0,±1 



An_ = 0,±I An z = or , . 
Au_ =0 An z = [ S) 



in a + and <r_ polarizations and: 

An+ = An- = An z = ±1 



(29) 



in 7T polarization. 

This means that a given n z ;fif) dressed eigen- 

state is connected to at most 2 states in the same Af 
ladder and 5 states in each of the Af ± 1 neighboring 
ladders. Among these 12 possible transitions, only those 
ending in a final state with lower energy are allowed in 
a spontaneous emission process. Using Eq. (^l|), we see 
that the energy loss is ±w z for the 2 7r-polarized tran- 
sitions, f2, J7 ± uj + , ± lj- for the 5 er+ transitions and 
— J7, — SI ± — Q ± uj_ for the 5 er_ transitions. From 
Eq. (p"5|), we know that lu + and oj~ are both smaller than 
n. Hence, all 5 transitions in <r_ polarization are forbid- 
den (by energy conservation) for spontaneous emission. 
On the other hand, all 5 transitions in er + polarization are 
allowed. In n polarization, only the transition An z = — 1 
is allowed by energy conservation. To summarize, in 
the harmonic approximation, a given |n+, n_, n z ) atomic 
eigenstate can only emit spontaneous photons at the six 
frequencies lo z , Q, D,±lu + , S1±oj_ the first transition being 
7T polarized, the 5 other ones having <r+ polarization. 

In the specific case where the initial state in the ground 
state wave packet |0, 0,0) with minimum values of the 
n +1 ri-,n z quantum numbers and maximum localization 
around the equilibrium point, all transitions decreasing 
any of the n+, n_, n z are obviously non-existent. Hence, 
we are left with only 3 possible transitions, all of them 
being cr+ polarized: 
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• Transition from the ground state wave packet 
|0,0,0;7V) to its image |0,0,0;7V- 1) at mi- 
crowave frequency SI. While the atomic state is 
not changed, a photon of the driving microwave is 
clastically scattered into another mode at the same 
frequency, but with a different direction. In the 
language of quantum optics, this process can be 
viewed as the elastic component in the resonance 
fluorescence of the driven atom. 

• Transition from |0, 0, 0; A/") to a n + — 1 dressed 
state 1 1, 0, 0; A/" — 1) occurring at frequency SI — uj+. 
This is an inelastic component. 

• Transition from |0,0,0;A/") to a n_ = 1 dressed 
state |0, 1,0; A/"— 1) occurring at frequency Sl + u_. 
This is also an inelastic component. Note that the 
difference on signs in the frequencies of the two in- 



elastic components, f2 — lo + and SI 



is directly 



related to the unusual - sign in the Hamiltonian, 
Eq. ©. 

In the harmonic approximation, we can go beyond the 
above qualitative predictions and calculate the partial 
decay rates to the different final states analytically. To 
this end, we need the matrix elements of the dipole op- 
erators between the initial and final states. This is easily 
obtained using the expansion of the dipole operator on 
the creation/annihilation operators in the normal modes, 
Eq. (p?j), as the matrix elements of these operators be- 
tween eigenstates are well known. We obtain: 



(0, 0, 0|ac - iy\0, 0, 0) = ST 2/3 <r 1/3 

-3q + Q- {q-Q)cn^ 

(i,o,o|^|o,o,o)= ' L>lJ '" 

(0,l,0|x-iy|0,0,0) = 



WQQN+ 
3q + Q- (q + Q)a^) 



(30) 
(31) 

(32) 



This finally makes it possible to calculate the decay 
rates of the 3 spontaneous emission processes, using 
Eq. @: 



r(0) = 



(33) 



for the elastic component, 
and 



|(l,0,0|a- ?2 /|0,0,0)| 2 (34) 



wJ) = 2q3( ° o + W ~ ? 1 (0, 1, 0\x - iy\0, 0, 0)i 2 (35) 



for the two inelastic components. The various quantities 
in these formula are related to the microwave frequency 
SI and to the scaled parameter q (itself related to the 
microwave amplitude by Eqs. ( |l0| ) and (0)) by Eqs. jl5| ) 



and ([32]). For a fixed classical dynamics (fixed q i.e. fixed 
scaled microwave amplitude Fq), the elastic decay rate 



scales as SI 5 / 3 



as does the decay rate of a circular 



state while the inelastic decay rates scale differently as 



si 2 = 



The results can be converted to usual units 



by multiplying them by me^/h 3 ~ 4.13 x 10 16 s -1 . 

In Fig. 0, we show the frequencies, squared dipole ma- 
trix elements and decay rates (in linear and logarith- 
mic scales) of the three components as a function of 
the scaled microwave field for principal quantum number 
no = 60, the corresponding microwave frequency being 
n/2% = 30.46 GHz. The first thing to note is obviously 
the dominance of the elastic component at the microwave 
frequency. About 80% of the spontaneous photons emit- 
ted by the system are at the microwave frequency. Elastic 
scattering of the microwave photons to other modes (dif- 
ferent direction) is by far the dominant process. The ra- 
diation diagram of these photons has the standard dipo- 
lar shape (this is identical for the three components as 
they are all due to the x — iy dipole operator). This 
dominance is easily understood from the different scal- 
ings of the three decay rates, Eqs. (|3|)-(|5|) discussed 
above. The dipole matrix element (0, 0, 0\x — iy\0, 0,0) 
measures the static dipole, i.e. the position of the cen- 
ter of the non-dispersive wave packet. It is thus - in the 
harmonic approximation - equal to x e , the coordinate of 
the classical fixed point, of the order of rig. In contrast, 
the non diagonal matrix elements (0,1, Q\x — iy|0,0,0) 
or (1,0, 0\x — iy\0, 0,0) measure how the excited states 
explore the vicinity of the fixed point. They are of the 

3/2 

order of n Q , i.e. about y/no times smaller. This is why 
the corresponding decay rates are much smaller. Because 
of the different scaling properties of the decay rates, it is 
easy to check that a change in the value of no gives sim- 
ilar plots for the decay rates of the inelastic components 
as a function of the scaled microwave field (they are just 
multiplied by a uniform numerical factor) while the elas- 
tic component is relatively increased in proportion with 
no. In other words, in the scmiclassical limit no — > oo, 
the inelastic components asymptotically vanish. 

Another important observation is that, although the 
two inelastic components have dipole matrix elements of 
comparable magnitude (see Fig. ||b), the decay rate of 
the SI — uj + component is several orders of magnitude 
smaller than the decay rate of the Sl+w^ component, see 
Fig. ||c. This is entirely due to the w 3 factor representing 
the density of electromagnetic mode. In practice, the 
SI — lo + component will be very difficult to detect. 

Finally, when the microwave field is sufficiently small, 
the various quantities can be expanded in powers of the 
scaled field. In this limit (see also Fig. ||d), the SI — uj + de- 
cay rate goes to zero as Fq whereas the Sl+w_ decay rate 
increases like 1/s/Fq and the elastic component SI keeps 
a constant decay rate. Of course, such a divergence in 
clearly unphysical and simply indicates the breakdown 
of the harmonic approximation. At lower field, the pen- 
dulum approximation of the Hamiltonian plfl has to be 
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used, restoring finite decay rates, roughly equally shared 
by the + c<j_ and Q components. 



iD„r 



2.0 | 

1.5 - 

1.0 — 

0.5 - 

0.0 — 

1.0 — 

0.5 r 

0.0 — 



£2-m 



(b) 



r[s" 



15 - 

10 - 

5 -* 

— 



(c) 



rof 1 ] io° 



(d) 



0.05 



0.10 



FIG. 2. Analytic predictions for the spectrum of spon- 
taneous photons emitted by a non-dispersive wave packet 
|0, 0, 0) of the hydrogen atom in a cicularly polarized mi- 
crowave of frequency Q./2n — 30.46 GHz corresponding to 
a principal quantum number no = 60, as a function of the 
scaled microwave held amplitude. Solid line corresponds to 
elastic scattering of the microwave photons, the dashed and 
dotted lines correspond to inelastic processes, (a): Frequen- 
cies of the emitted photons, (b) : Squares of the dipole matrix 
elements involved in the 3 transitions. The elastic component 
is largely dominant, (c): Decay rates along the 3 transitions. 
Because of the rapid increase in the density of electromag- 
netic modes with frequency, the fi — ui+ component is almost 
invisible, (d): Same as (c) on a logarithmic scale. The decay 
rates of the — ui+ , Q and £1 + uj- components roughly scale 



as Fq n , n and F 



-V2„ 



respectively. 



IV. COMPARISON WITH EARLIER RESULTS 

The former approach to spontaneous emission of the 
non-dispersive wave packets has been carried out 
in the rotating frame using the Hamiltonian Eq.(^), or 
rather its generalization which includes also other non 



occupied electromagnetic modes. Using the notion of ro- 
tational frequency shift found previously by the same au- 
thors , the frequency of the emission in the laboratory 
frame has been found. A single frequency, since the au- 
thors considered the emission at Q + lj_ only assuming, 
incorrectly, that spontaneous emission of a photon in the 
rotating frame must lead to a decrease of the energy (as 
it happens in the laboratory frame). 

To see the fallacy of this argument, consider the stan- 
dard spontaneous decay of the circular |n+l, n, n) hydro- 
genic state in the absence of any driving electromagnetic 
field. In the dipole approximation, it can decay sponta- 
neously only to the |n, n— 1, n— 1) circular state. Consider 
the same situation in the frame rotating at frequency f2. 
The Hamiltonian in the rotating frame is obtained from 
Eq. (|) with F = 0. The energy of the initial state in the 
rotating frame is then Ei = — 1/2 (n+ l) 2 — Vln while the 
energy of the final state Ef = — l/2n 2 — £l(n— 1) is greater 
than Ei provided > (n + l/2)/(n(n + l)) 2 « 1/n 3 . 
Hence, emission of the spontaneous er + polarized photon 
increases the energy in the rotating frame. 

The same situation occurs for the standard toy model 
of quantum optics: the two-level atom interacting with 
a quasi-resonant monochromatic field, treated in the ro- 
tating wave approximation. In the "rotating frame" p3| , 
spontaneous emission does not occur only from the up- 
per to the lower state of the (Autler-Townes) doublet. 
This particular line (upper to lower state in the doublet) 
occurs at the frequency of the driving field plus the Rabi 
frequency. But the emission may also take place from 
the lower to the upper state of the doublet (at driving 
frequency minus the Rabi frequency) as well as from the 
lower (upper) state to "itself . In effect the resonance flu- 
orescence spectrum has the famous Mollow triplet struc- 
ture p2|p| with the dominant elastic component in the 
middle. 

This discussion shows that the rotating frame approach 
leads to quite counter intuitive picture. By comparison, 
the dressed state picture presented in the previous sec- 
tion seems to be quite intuitive. One may easily verify 
that extending the calculation presented in to include 
the lines emitted at the microwave frequency f2 (which 
corresponds in the rotating frame to the emission from 
the wave packet state to itself) as well as at fi — w+ leads 
to the same results as those obtained in the dressed atom 
picture above |33| . This shows that both approaches lead, 
as it should be, to the same spontaneous emission rates. 

From a practical point of view, as the decay at £1 — lo + 
is much slower than the decay at + lj- (see discussion 
above), the numbers given in can be considered correct 
f{4| if elastic scattering does not have to be taken into 
account. A non-trivial question is whether elastic scat- 
tering of photons has to considered as a decay process or 
not. In our opinion, this is simply a question not properly 
formulated. Thinking resonance fluorescence in terms of 
decay is probably inadequate as the system is continu- 
ously fed with external photons. If one measures the 
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population in the ground state wave packet |0, 0, 0), this 
quantity will certainly decay because of inelastic com- 
ponents only. If one measures the spectrum radiated by 
the atom, the answer is completely different as the elastic 
component dominates. Also, if one considers the widths 
of the various components, the elastic process will come 
into the game as it destroys the phase coherence of the 
initial wave packet and consequently broadens the line. 
A general discussion of such effects and more generally 
of resonance fluorescence in the framework of quantum 
optics can be found in p2fl . 

Let us also note that the "classical" emission rate of 
an electron moving on a circular orbit presented in Q| is 
incorrect. A classical charge moving on a circular orbit 
with frequency Q must emit photons at the same fre- 
quency. In our opinion, the only meaningful quantity is 
the decay rate of the energy in the atomic system. Clas- 
sically, it is proportional to the square of the acceleration 
and, for the classical periodic orbit at the microwave fre- 
quency, it is in atomic units |p5[: 



dE 
~dt 



2a 3 ft 8 /V 2/3 



(36) 



Considering that, in the semiclassical limit % -> oo, the 
system radiates essentially photons as frequency fi, we 
immediately see that this classical rate exactly (i.e. with- 
out additional numerical factor of the order of one) co- 
incides with the quantum decay rate of energy deduced 
from Eq. (33). It also exactly corresponds to the deca' 
rate of atomic circular states, as given in e.g. Ref. J56 
in the limit of vanishing microwave field q — > 1. In 
the decay rate in the Q + cj_ component is incorrectly 
compared to the relative decrease of the classical energy 
in time. The first quantity - as discussed above - is 
roughly no times smaller than the elastic decay rate while 
the relative decrease of classical energy in time is also 
about hq times smaller than the rate of emission of pho- 
tons because the energy of a radiated photon is about no 
times smaller than the total energy. The two no depen- 
dences cancel out, but not the numerical factors in front. 
This explains why a false classical-quantum correspon- 
dence (within constant numerical factors) is observed in 
. Both classical and quantum decay rates are - in our 
opinion - underestimated by a ~ no factor. 



HOW GOOD IS THE HARMONIC 
APPROXIMATION ? 



From the beginning of investigations on non-dispersive 
wave packets in the hydrogen atom driven by a circu- 
larly polarized microwave, there has been a controversy 
regarding the accuracy of the harmonic approximation 
used when introducing the wavepackets as Gaussians cen- 
tered around the equilibrium point |l]], see Section II. 
Firstly in a Comment |l5[ ] and later in a number of pa- 
pers ]l6|-pcfl Farrelly and coworkers claimed that "highly 



anharmonic nature of the effective potential in the vicin- 
ity of the equilibrium point will cause any such putative 
coherent state to disperse .." unless the wave packet is 
built from states of unrealistically high quantum num- 
ber (no > 200). They have shown that an additional 
static magnetic field applied in the direction of microwave 
propagation (Oz axis) allows to relatively decrease the 
importance of the anharmonic terms in the expansion 
around the equilibrium point. Further they claim that 
"because these extrema are locally harmonic their vac- 
uum states are truly coherent states in the original sense 
of Schrodinger" [g(J (by extrema one should understand 
global equilibrium points). 

Our understanding of the phenomenon is slightly dif- 
ferent. While Farrelly and coworkers stress that the 
Gaussian shape of the wave packet is essential for non 
spreading time evolution in the locally harmonic poten- 
tial, we adopt a bit less narrow-minded point of view. 
We have shown jl| that, even if the harmonic approxi- 
mation is only partially valid, there will be a one-to-one 
correspondence between the quasienergy states of the ap- 
proximate harmonic Hamiltonian, Eq. (vw, and the full 
Hamiltonian in the rotating frame, Eq. (H). In particu- 
lar, the ground state of the local harmonic approximation 
(referred to above as ground state wave packet) is a good 
approximation to an exact eigenstate of the full problem 
(or via the unitary transformation to the corresponding 
Floquet or dressed state). We refer to this exact eigen- 
state of the full atom+microwave field problem as the 
non-dispersive wave packet since: 

• it is a linear combination of atomic states; 

• it is centered at the stable equilibrium point in the 
rotating frame, thus it follows the classical circu- 
lar periodic orbit of the problem in the laboratory 
frame; 

thus fulfilling standard textbook criteria for "a wave 
packet" . Being an eigenstate of if ro t , Eq. (B) , the state 
is rigorously non spreading provided we neglect the cou- 
pling to the continuum - i.e. the possibility of ionization. 
We have shown ||-[3| that typical lifetimes for the ioniza- 
tion exceed many thousands of Kepler periods (reaching 
up to 10 6 Kepler periods for no ~ 60 or more). 

In our understanding of the phenomenon, the validity 
of the harmonic approximation is not crucial for the very 
existence of non-dispersive wave packets in the Hydrogen 
atom driven by a circularly polarized microwave. More- 
over, the phenomenon is more general and, as mentioned 
in the introduction, may be found in other systems such 
as atoms driven by a linearly polarized microwave 
or one-dimensional models |25| supporting stable classi- 
cal 1:1 resonance islands. In the general case, the wave 
packet preserves its shape at intervals equal to the period 
of the external driving, only in the circular polarization 
problem (due to the additional symmetry which allows 
for the complete removal of the time-dependence) the 
wave packet shape is the same at all times. There is no 
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reason to restrict the notion of non-dispersive (or non 
spreading) wave packet to a single particular situation 
which can be made almost fully harmonic. One should 
also remember that even in the presence of a magnetic 
field, a case favored by Farrelly and coworkers, the poten- 
tial is also only approximately harmonic, and Gaussian 
wavepackets will eventually disperse. While wave packets 
in our sense, being single dressed states, will not. In any 
case, the addition of a magnetic field does not change the 
nature of the problem, it may just keep some anharmonic 
terms smaller. 

While not necessary for the very existence of non- 
dispersive wave packets, the harmonic approximation is 
useful for analytic evaluation of different properties, such 
as spontaneous emission discussed above. In view of Far- 
relly and coworkers remarks, it would seem that the re- 
sults obtained in the previous section are incorrect. The 
simplest way to check this statement is to evaluate the 
"exact" spontaneous emission spectrum from the ground 
quasicnergy cigcnstatc of the full Hamiltonian. 

The presence of the anharmonic terms in the full 
Hamiltonian has to manifest itself in the break-up of the 
selection rules discussed in Section II. Instead of three 
possible lines in the o~ + emission, one should expect a 
multitude of lines. Using the harmonic level assignment, 
the |0, 0, 0; N) wave packet may then decay to e.g. states 
approximated by k;fif — 1) with k being larger 
than unity (with k even since o~ + dipole element may link 
states of same z-parity only). Also, due to couplings in- 
troduced by nonlinear terms, spontaneous emission with 
other (i.e. cr_ or tt) polarizations becomes possible. 

Thus, the spontaneous emission spectrum provides a 
useful test of the validity of the harmonic approximation. 
Having at our disposal the computer code which allows 
reliably to find the exact quasienergy states of the prob- 
lem ^-|(| , it has been sufficient to evaluate dipole matrix 
elements between the exact quasienergy states. The re- 
sulting squared dipole moments are shown in Fig. |^ to- 
gether with the corresponding squared dipoles as given 
by the harmonic approximation, Eq. (|32|), presented as 
three crosses. Observe that, the harmonic approximation 
works quite well. The selection rules are well preserved, 
the frequencies of the emission (line positions) agree very 
well, the dipole values (heights of the lines) agree within 
10% for all three lines. This shows that the harmonic 
approximation is quite good, in the absence of any mag- 
netic field, not only for the ground state wave packet, 
|0,0,0) but also for primary lo + [state |1,0,0)] and w_ 
[state |0, 1, 0)] excitations. Observe that the height of the 
peak at fl is proportional to |(0, 0, 0|x — iy\0, 0,0)| 2 in- 
volving the wavefunction of the ground state wave packet 
only. The sidebands, corresponding to transitions to ex- 
cited wave packets provide independent tests of the cor- 
responding wavefunctions. 

Fig. |] shows the corresponding spontaneous emission 
rates. The component at — u> + is not visible due to the 
low frequency of the emission. On the other hand, the 
spectrum at high frequencies is enhanced by cubic power 



of the frequency factor coming from the photonic density 
of states. Here, careful inspection of the spectrum cor- 
responding to exact eigenstates reveals a tiny peak cen- 
tered at Q + 2uj- (indicated by an arrow in the figure). 
Such a transition is forbidden in the harmonic approxi- 
mation, its presence indicates that the harmonic approx- 
imation is not exact. Still, the smallness of anharmonic 
terms strongly indicates that the harmonic approxima- 
tion works well, at least for the typical value no = 60. 

If one considers lower no (say around no = 30) when 
the effective ft, is bigger, the anharmonic terms in the 
Hamiltonian become more important and the sponta- 
neous emission spectrum contains other lines, at £l + ku)- 
with k > 2 or combinations, e.g. at Q + 2lj- — lu + . This 
nicely correlates with our study of ionization mechanism 
||. Starting around no = 40, ionization clearly occurs 
via a chaos assisted tunneling mechanism (obeying ex- 
ponential behavior with ft change) with clear deviations 
from such a behavior for lower no. 

Similarly, much weaker microwave amplitudes (say 
Fq = 0.005) when the resonance island is so small that 
the quantum eigenstate does not fit into it, lead to a sim- 
ilar behavior of the spontaneous emission spectrum - the 
presence of overtones. Still we have checked that in quite 
a broad range of 0.01 < F < 0.06, the harmonic ap- 
proximation works similarly to the case shown in Fig. ^ 
and Fig. |], provided no accidental avoided crossing with 
other Floquet states occurs |37j. 

While the calculations presented in this paper are con- 
sistently done at the lowest order of the harmonic approx- 
imation using Cartesian coordinates, similar calculations 
can be done using, for example, the spherical coordinates 
pT| , p2| . Consistently, the results at lowest order are of 
course identical. However, it is known that - within the 
harmonic approximation - the wavefunction in spheri- 
cal coordinates is a better approximation to the "exact" 
wavefunction than the corresponding one in Cartesian 
coordinates, as it incorporates some of the anharmonic 
terms; the price to pay is the breaking of the selection 
rules, Eqs. (p8|)-(]2^). Hence, the transition rates calcu- 
lated using harmonic wavefunctions in spherical coordi- 
nates incorporate already some of the higher order terms 
and might be closer to the "exact" results. Note, how- 
ever that if the calculation is performed consitently at 
any order in perturbation theory (the anharmonic terms 
being the perturbation), the results should of course be 
independent of the choice of coordinates. 
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FIG. 3. Single photon spontaneous emission from the 
non-dispersive wave packet of the hydrogen atom in a cicu- 
larly polarized microwave of frequency Q = 1/60 3 (in atomic 
units) corresponding to a principal quantum number no = 60, 
and scaled amplitude Fq — 0.04445881 (corresponding to 
q = 0.9562). We have here represented the square of the 
scaled dipole matrix element Do = D jn\ connecting the ini- 
tial |0, 0, 0) wave packet to the final state as a function of 
the scaled energy difference (hence, ratio of the frequency of 
the emitted photon to the microwave frequency fi) between 
the two states. The stick spectrum is the "exact" result as 
obtained from a numerical diagonalization and the crosses 
represent the analytic result obtained using the harmonic ap- 
proximation, Eq. p2]). There are three dominant lines at fre- 
quencies Q — UJ+, fl and fi + w_, all of them a + polarized. The 
strongest line (by a factor proportional to no) is the elastic 
scattering of the microwave photons by the dressed atom. The 
other transitions (as well as transitions with a_ or n polariza- 
tions) are negligible, the corresponding sticks have practically 
zero heights (yielding a broken line build of individual sticks 
each representing one of the over 400 Floquet final states), 
which proves the validity of the harmonic approximation. 



FIG. 4. Same as figure |3|, but for the decay rates along the 
different transitions and in S.I. units. The density of modes of 
the electromagnetic wave introduces an additional u 3 factor 
which completely kills the transition at frequency fl — 
invisible in the figure. One can see a small line at frequency 
approximately £1 + 2lj_ (see arrow), an indication of a weak 
breakdown of the harmonic approximation. 

VI. CONCLUSIONS 

It is worth stressing again that, in the semiclassical 
limit (small Q or large no), the dominant contribution 
to the spontaneous emission comes from the elastic com- 
ponent at frequency Q which does not destroy the wave 
packet but rather leads to the decoherence in its time- 
evolution. The wave packet converts the energy of the 
driving photon into a spontaneously emitted photon, 
moving at the same time on the circular orbit. 

From a practical point of view, the spontaneous emis- 
sion rates are several orders of magnitude smaller than 
the ionization rate in a typical situation no = 60 and 
Fq = 0.05 0, and thus unlikely to be observable. How- 
ever, this ratio can be made smaller by increasing no- 
Indeed, as shown in |(|, the typical ionization rate de- 
creases exponentially with no (because it is due a chaos 
assisted tunneling process) while the spontaneous decay 
rates decrease algebraically with no. For very large no, 
spontaneous emission will eventually become the domi- 
nant process. For Fq = 0.05, the cross-over should take 
place around no = 200 which makes it very difficult to 
observe experimentally. Another interesting possibility 
is to decrease the scaled field Fo so that the ionization 
rate drops considerably. According to numbers shown in 
, the cross-over for no = 60 should take place around 
Fq = 0.03 in a regime where the harmonic approxima- 
tion is still valid (the results are not shown since they 
are very similar to those depicted in Fig. ^). In any case, 
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the total decay rate is of the order of few tens of photons 
per second, making the experimental observation rather 
difficult. 
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